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HIGHLIGHTS 


• Observations of ozone are higher over the Chesapeake Bay than areas upwind on land. 

• Dry deposition rates, boundary layer depth, and photolysis play an integral role. 

• Model resolution plays a role in determining accurate surface ozone concentrations. 

• Observations of total reactive nitrogen are much lower than model simulations. 
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Air quality models, such as the Community Multiscale Air Quality (CMAQ) model, indicate decidedly 
higher ozone near the surface of large interior water bodies, such as the Great Lakes and Chesapeake Bay. 
In order to test the validity of the model output, we performed surface measurements of ozone (O3) and 
total reactive nitrogen (NO y ) on the 26-m Delaware II NOAA Small Research Vessel experimental (SRVx), 
deployed in the Chesapeake Bay for 10 daytime cruises in July 2011 as part of NASA’s GEO-CAPE CBODAQ 
oceanographic field campaign in conjunction with NASA’s DISCOVER-AQ air quality field campaign. 
During this 10-day period, the EPA 0 3 regulatory standard of 75 ppbv averaged over an 8-h period was 
exceeded four times over water while ground stations in the area only exceeded the standard at most 
twice. This suggests that on days when the Baltimore/Washington region is in compliance with the EPA 
standard, air quality over the Chesapeake Bay might exceed the EPA standard. Ozone observations over 
the bay during the afternoon were consistently 10—20% higher than the closest upwind ground sites 
during the 10-day campaign; this pattern persisted during good and poor air quality days. A lower 
boundary layer, reduced cloud cover, slower dry deposition rates, and other lesser mechanisms, 
contribute to the local maximum of ozone over the Chesapeake Bay. Observations from this campaign 
were compared to a CMAQ simulation at 1.33 km resolution. The model is able to predict the regional 
maximum of ozone over the Chesapeake Bay accurately, but NO y concentrations are significantly over- 
estimated. Explanations for the overestimation of NOy in the model simulations are also explored. 

© 2013 Elsevier Ltd. All rights reserved. 


1. Introduction 

Surface ozone, in high enough concentrations, is a hazardous 
secondary air pollutant regulated by the United States Environ- 
mental Protection Agency (U.S. EPA). As of 2013, the U.S. EPA has 
set an ambient 8-h daily maximum concentration of 75 parts per 
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billion by volume (ppbv) under the National Ambient Air Quality 
Standards (NAAQS). This standard has been set based on several 
health studies (Anenberg et al., 2010; Bell et al., 2004; Fann et al., 
2011) that show inhaling ozone can lead to premature aging of the 
lungs, difficulty breathing, increased risk of asthma attacks, and in 
rare cases death. Reactive nitrogen species and volatile organic 
compounds, emitted by natural and anthropogenic sources, react 
photochemically to create ozone (Seinfeld and Pandis, 2006). 
Warm temperatures, along with a capping inversion, associated 
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with a strong summer anticyclone aid in the rapid growth of 
ozone near the surface during the late morning and early 
afternoon. 

Air quality models such as the Community Multiscale Air 
Quality (CMAQ) model indicate decidedly higher ozone near the 
surface of large interior waters bodies such as the Great Lakes and 
Chesapeake Bay (e. g., Godowitch et al., 2008). In order to test the 
validity of the model output, we performed surface measurements 
of ozone (O 3 ) and total reactive nitrogen (NO y ) on 26-m Delaware II 
NOAA Small Research Vessel experimental (SRVx), deployed in the 
Chesapeake Bay for ten daytime cruises in July 2011. The objectives 
of this paper are to: 

• Compare ozone observations over the Bay to nearby land areas 

• Determine if ozone concentrations are indeed higher over the 
Bay 

• Determine if known meteorological and chemical processes can 
explain the observed differences 

• Investigate whether model grid resolution plays a role in 
determining the simulated surface ozone concentrations over 
the Bay 

• Investigate NO y observations to determine if this group of pre- 
cursors is accurately predicted by the model simulations 


1.1. Current EPA regulations and designations 

The U.S. EPA has designated a “moderate” non-attainment re- 
gion in the counties surrounding Baltimore, MD and a “marginal” 
non-attainment region throughout the New York City- 
Philadelphia-Washington, D.C. metropolitan corridor (U.S. EPA, 
2008). A “moderate” attainment region is an area that has a 
ground station that exceeds an 8 -h maximum concentration of 
85 ppbv O3 in its 3rd highest day, while a marginal attainment 
region has a threshold of 75 ppbv. 

The Baltimore /Washington metropolitan region exceeds the 
threshold due to favorable meteorological conditions for ozone 
production during the summer months and substantial ozone 
precursor emissions, generated locally as well as advected to the 
region during strong westerly transport conditions (Ryan et al., 
1998; He et al., 2013). Peaks in surface ozone are highest just 
downwind of major metropolitan areas due to the enhanced 
emissions from the metropolitan city centers (Kleinman et al., 
2000). This has been shown in many air quality model simula- 
tions (Yegorova et al., 2011; Castellanos et al., 2011) and has been 
verified by ground monitoring stations (U.S. EPA, 2006; Castellanos 
et al., 2009). In the Baltimore-Washington region there are com- 
plex interactions that arise with the influence of the Chesapeake 
Bay breeze (Loughner et al., 2011 ; Stauffer et al., 2012), which have 
not been fully investigated. 


1.2. Previous field campaigns over interior water bodies 

Measurements of ozone over the southern Great Lakes during a 
2007 summer field campaign show higher concentrations of ozone 
over the lakes than over the adjacent land with the biggest differ- 
ence detected at night (Levy et al., 2010). A similar study was 
conducted over Lake Michigan in the summers of 1990 and 1991, 
where O3 and NO x were monitored from aircraft (Luria et al., 1992). 
High levels of ozone are shown only at the lowest levels of the 
boundary layer, which they attribute to a lack of vertical mixing 
over the lake (Dye et al., 1995). An experiment in 2003 measured 
ozone at the Chesapeake Bay Lighthouse, located on an island 15 
miles to the east of the entrance to the Chesapeake Bay, as a means 
to test ozone monitoring on ocean buoys and towers (Hintsa et al., 
2004). This field campaign found ozone at the surface consistently 
exceeding 80 ppbv during an air quality episode from June 24 to 28, 
2003. 

1.3. DISCOVER- A(l And GEO-CAPE CBODAQ field campaigns 

During the month of July 2011, the National Aeronautics and 
Space Administration (NASA) conducted a comprehensive air 
quality field study, DISCOVER- AQ( Crawford and Pickering, 2011 ), in 
the Washington, DC-Baltimore, MD metropolitan area and over the 
Chesapeake Bay to investigate air quality with the primary goal of 
providing data to better interpret observations from current and 
future satellites for air quality applications. In conjunction with 
DISCOVER-AQ NASA conducted the oceanographic field campaign 
GEO-CAPE CBODAQ (Geostationary Coastal and Air Pollution 
Events-Chesapeake Bay Oceanographic Campaign with DISCOVER- 
AQ), to address questions related to both estuarine biogeochemical 
processes as well as atmospheric pollution over the Chesapeake 
Bay urban estuarine environment (Tzortziou et al., 2013). There 
were seven instrument platforms during the field campaign. A 
detailed description of the modes, locations, types, and days of 
observations is provided in Table 1. This paper focuses on obser- 
vations from the 26 m Delaware II NOAA Small Research Vessel 
experimental (SRVx) deployed in the Chesapeake Bay as part of the 
CBODAQ campaign from July 11 to 20, 2011. 

2. Materials and methods 

2.1. Measurements description 

The SRVx was equipped with a Thermo Environmental Model 49 
UV photometric ozone (O3) analyzer and a modified Thermo 
Environmental Model 42C chemiluminescence nitric oxide (NO) 
analyzer retrofitted with an external molybdenum catalyst to also 
measure total reactive nitrogen (NO y ) (Delany et al., 1982). The NO y 
analyzer was zeroed for 10 min each hour during the campaign and 


Table 1 

Modes of measurement during Phase I of the DISCOVER-AQ. campaign in the Baltimore-Washington region. 


Mode of measurement 

Trace gases measured 

Location 

Days in which Active during 
July 11, 2011 through July 20, 2011 

NASA P3-B 

0 3) NO, N0 2 , PNs, ANs, HN0 3 , 
NOy, CO, CH 4 , CH 2 0, C0 2 , VOCs 


July 11, July 14, July 16, July 20 

UC-12 king air 

Column 0 3 , N0 2 , CH 2 0 


July 11, July 14, July 20 

Cessna 402B 

0 3 , S0 2 , CO, no 2 


July 11, July 18, July 20 

MDE ground stations 

0 3 , NO, NOy, PM2.5, CO, S0 2 

Aldino, MD, Beltsville, MD, Edgewood, MD, 
Essex, MD, Fairhill, MD, Padonia, MD 

Every day 

Pandora 

Column N0 2 

Aldino, MD, Beltsville, MD, Edgewood, MD, 

Essex, MD, Fairhill, MD, Padonia, MD, 

College Park, MD, Catonsville, MD, SRVx NOAA Vessel 

Every day 

Ozonesondes 

O 3 

Edgewood, MD 

Every day 

SRVx NOAA research vessel 

0 3 , NO, total NOy 

See Fig. 1 

Every day 
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measurements were adjusted based on the drift of the instrument. 
The NO y analyzer was calibrated in-situ on July 19, 2011 using a NO 2 
standard reference material (SRM) from the National Institute of 
Standards and Technology (NIST). 

2.2. Model description 

In this study, we use U.S. EPA’s Community Multiscale Air 
Quality (CMAQ) (Byun and Schere, 2006) model Version 5.0, driven 
off-line by output from the Weather Research and Forecasting 
(WRF) (Skamarock et al., 2008) model Version 3.3 to simulate the 
state of the atmosphere covering the entire months of June and July 
2011. Passing the meteorology into an air quality model at a high 
temporal resolution or running the chemistry online within a 
meteorological model is preferable, but requires significantly more 
computational resources (Grell et al., 2004). 

The WRF and CMAQ model simulations are at 36, 12, 4 & 1.33 km 
resolution in the area of interest with 34 verticals levels from the 
surface to 100 mbar and 16 levels within the lowest 2 km in order to 
accurately simulate boundary layer processes. The 1.33 km model 
domain covers the Baltimore— Washington metropolitan region 
and nearby Chesapeake Bay. The North American Regional Rean- 
alysis (NARR) is used for the model initial and outermost lateral 
boundary conditions in WRF. The Multi-scale Ultra-High Resolution 
(MUR) dataset was used to set the sea surface temperatures. The 
WRF model was re-initialized every 3 days and run in 3.5 day in- 
crements. The first 12 h of each simulation was thrown out (i.e., not 
passed to CMAQ). WRF model output is input into the Meteo- 
rology-Chemistry Interface Processor (MCIP; Otte and Pleim, 2010) 
to create meteorological input fields for CMAQ. Chemical initial and 
boundary conditions come from a MOZART-4 simulation (Emmons 
et al., 2010). The Carbon-Bond-05 (CB05) gas-phase chemical 
mechanism (Yarwood et al., 2005) was used in CMAQ. The CMAQ 
and WRF simulations began May 24, 2011, which allows ample 
spin-up time for our comparison in mid-July. 

Anthropogenic emissions input files for CMAQ are created with 
the Sparse Matrix Operator Kernel Emissions (SMOKE) modeling 
system (Houyoux and Vukovich, 1999). We use a projected 2012 
emissions inventory because a 2011 emissions inventory is not yet 
available. Annual projected point and countywide area emissions 
are temporally distributed based on the time of day, day of the 
week, and season based on temporal surrogates from the EPA. 
Mobile emissions estimates from cars, trucks, and motorcycles are 
computed with the Motor Vehicle Emission Simulator (MOVES; 
Kota et al., 2012). Point sources are vertically distributed based on 
the meteorology, stack height, and the temperature and velocity of 
the emissions exiting the stack. Biogenic emissions are calculated 
using Biogenic Emissions Inventory System (BEIS) and lightning 
NO x emissions are calculated in-line within the CMAQ model. 
Further model details on the parameterization options are 
described in the Supplementary Material. 

3. Results 

3.1 Observational comparisons: ozone 

The SRVx was deployed in the Chesapeake Bay for 10 daytime 
cruises during the DISCOVER-AQ Maryland campaign extending 
from July 11, 2011 through July 20, 2011. This overlapped with four 
flights of the NASA P3-B (a four-engine turboprop capable of long 
duration flights of 8—12 h) three flight days (2 flights per day) of the 
UC-12B King Air (a twin-engine turboprop capable of 6 h flights) 
and three flight days (2 flights per day) of the University of Mary- 
land (UMD) Cessna 402B (a twin-piston engine, unpressurized 
aircraft) (Table 1 ). The SRVx docked each night in Annapolis, MD 


and had different cruise route each day (Fig. 1). The instruments 
were running while the SRVx was in port overnight in Annapolis, 
MD, but the data are subject to frequent local emissions. 

A time series of O3 for the 10-day period can be seen in Fig. 2. On 
four days ozone exceeded the 8-h maximum 75 ppbv NAAQS 
threshold on the moving vessel in the Chesapeake Bay: July 12, 13, 
19 & 20. During this same time period, ground stations in the 
DISCOVER-AQ field campaign region (stations denoted in Table 2) 
exceeded the 75 ppbv threshold an average of 0.71 times per 
ground station. This alone is an indicator that the ozone may be 
higher near the surface of the Chesapeake Bay than nearby ground 
stations. 

Comparing the hourly ozone at the SRVx’s location and closest 
upwind ground station reinforces the idea that higher ozone con- 
centrations exist over the Bay. The closest upwind ground station 
was determined by using the backward trajectories at 10 m, 500 m, 
1500 m heights ending at 1800 UTC (2 PM local time) using Global 
Data Assimilation System (GDAS) meteorological data in the NOAA 
HYSPLIT trajectory model (Draxler and Rolph, 2003). Eight-hour 
maximum ozone from all relevant ground sites and the SRVx can 
be seen in Table 2. The closest upwind ground station was often 
only 20-30 km away and was chosen to ensure that comparisons 
were made following the same parcel of air. Fig. 3 a and b shows that 
during an exceedance day (July 13) and non-exceedance day (July 
14) on the boat, the ozone near the surface of the Chesapeake Bay is 
uniformly higher. During the afternoon of July 13 the ozone mea- 
surement on the SRVx was 10-40 ppbv greater than at the Calvert 
County MDE site. Ozone was consistently 10-20 ppbv greater over 
the Bay than at the Essex MDE site throughout the day on July 14. 

The 8-h maximum ozone concentration over the Bay during 
each day of the 10-day cruise averaged 12.7 ± 6.1 ppbv higher than 
the closest upwind ground site. The systematic high anomaly over 
the Chesapeake Bay can be seen in Fig. 4. The closest upwind 
ground site never experienced higher 8-h maximum ozone and 



Fig. 1. Map of NOAA Delaware II SRVx routes from July, 11, 2011 through July 20, 2011. 
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Date of July 2011, tickmark centered at 12 PM local time 

Fig. 2. Ozone concentration (ppbv) as a function of time from July 11, 2011 through 
July 20, 2011. Map routes for each specific day can be seen in Fig. 1. From 7 PM until 6 
AM local time, the boat was docked at the US Naval Academy in Annapolis, MD. 



Fig. 3. a. Ozone concentration on July 13, 2011 (ppbv) as a function of time at the 
SRVx’s location and the Calvert County ground monitoring station, the closest upwind 
monitoring station, b. Same as Fig. 3a but on July 14, 2011 Essex was the closest upwind 
monitoring station. 


only during three days did any ground station in the region have an 
8-h maximum ozone concentration that was 10 ppbv higher than 
the SRVx’s location. This was especially pronounced on July 13 
when the SRVx saw an 8-h maximum of 85 ppbv and none of the 
ground stations in the region exceeded the 75 ppbv NAAQS stan- 
dard. When compared to 8-h maximum ozone at the ground sta- 
tions in the Baltimore “moderate” non-attainment area, the 8- 
h ozone at the SRVx’s location was 4.6 ± 14.3 ppbv higher sug- 
gesting that the Chesapeake Bay has just as poor if not worse air 
quality than the surrounding “moderate” non-attainment area. 

The ozone concentration remained higher over the Chesapeake 
Bay later into the afternoon than over the ground stations, sug- 
gesting that there must be a mechanism to maintain high O3 con- 
centrations later into the day. A plot of the median hourly ozone 
concentrations at the SRVx’s location and closest upwind ground 
station (Fig. 5) illustrates the late afternoon high anomaly. Ozone 
concentrations over the Bay are greater and exist for longer dura- 
tions than over the upwind land area due to several potential 
causes: 

(1) A difference in ozone deposition rates over land and water; 

(2) A shallower PBL depth over the Chesapeake Bay than the 
nearby land causing emissions from shipping to be trapped 
near the surface; 

(3) Fewer fair-weather cumulus clouds over the Chesapeake Bay 
allowing for increased photolysis; and 


(4) Decreased boundary layer venting caused by a meso-high 
pressure that develops over the Bay as part of the bay- 
breeze circulation trapping pollutants. 

Furthermore, when meteorological conditions are conducive, a 
low-level jet can form overnight transporting polluted air over the 
Chesapeake Bay from the Norfolk/Virginia Beach, VA metropolitan 
region bypassing ground stations allowing for increased ozone 
production over the Bay. This phenomenon, however, was not 
observed during this field campaign. 

3.2. Model comparisons: ozone 

Our CMAQ model simulation results typically reproduce the 
systematically higher ozone concentrations over the Chesapeake 
Bay than in the Baltimore-Washington region. As shown in Fig. 6, 
median ozone concentrations for the 10-day period output by both 
the 1.33 km and 4 km resolution CMAQ model simulations closely 
match the observations from the SRVx throughout the day. Model 
mean bias of ozone at the boat’s location was 0.78 ppbv, but a root- 
mean square error (RMSE) of 10.14 ppbv Table 3 shows the mean 
model bias, normalized mean bias, root-mean square error, and 
normalized mean error. At a grid cell size of 12 km, the surface 
ozone output by the model begins to lose correlation and at a grid 
cell size of 36 km, there was very little correlation throughout the 


Table 2 

Maximum 8 h ozone at various sites in the Maryland & Delaware region. 


Site name 

7/11/11 

7/12/11 

7/13/11 

7/14/11 

7/15/11 

7/16/11 

7/17/11 

7/18/11 

7/19/11 

7/20/11 

SRVx Boat 

58 

77 

85 

60 

50 

51 

52 

66 

80 

86 

Aldino 

69 

66 

56 

51 

70 

59 

61 

88 

73 

77 

Calvert Co. 

58 

71 

67 

49 

38 

37 

40 

57 

66 

59 

Davidsonville 

53 

67 

61 

50 

45 

41 

41 

81 

76 

63 

Edgewood 

61 

61 

52 

49 

66 

55 

53 

82 

59 

74 

Essex 

58 

61 

53 

51 

58 

47 

51 

71 

72 

68 

Fairhill 

68 

65 

49 

48 

69 

51 

55 

86 

68 

66 

HU-Beltsville 

75 

63 

59 

61 

52 

49 

54 

75 

60 

80 

Millington 

59 

68 

58 

54 

59 

46 

47 

62 

68 

68 

Padonia 

73 

64 

54 

52 

62 

60 

53 

67 

65 

86 

PG Equestrian 

54 

72 

62 

51 

44 

42 

41 

80 

75 

61 

South Maryland 

55 

70 

63 

46 

40 

38 

42 

54 

64 

71 

Lums Pond, DE 

59 

65 

53 

49 

58 

N/A 

49 

76 

61 

57 

Seaford, DE 

53 

76 

62 

44 

46 

N/A 

41 

52 

63 

65 

Closest ground site 

58 

67 

62 

49 

38 

37 

40 

54 

64 

68 

Name of closest ground site 

Calvert 

Davidsonville 

PG Equestrian 

Lums Pond Calvert 

Calvert 

Calvert 

S Maryland 

S Maryland 

Millington 
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Fig. 4. 8-h maximum ozone concentrations (ppbv) at the SRVx’s location and the 
closest upwind ground monitoring station from July 11, 2011 through July 20, 2011. 

day; both the 12 km and 36 km model runs show a high model bias 
in the late morning and afternoon. Model resolution seems to play 
an integral role in predicting ozone concentrations over the Bay. 

Furthermore, we conducted a validation of the 2-m temperature 
in the 1.33 km WRF model simulation to ensure that the meteo- 
rology is indeed representative of the actual conditions during the 
10-day period. Model mean bias of 2-m temperature at the boat’s 
location over the 10-day period was -0.52 °C and a RMSE of 1.59 °C, 
which was a lower error than the nearby BWI airport. Table 4 shows 
the hourly mean model bias and RMSE at the boat’s location and 
the BWI airport. WRF at 1.33 and 4 km resolution was able to 
reasonably capture wind speed and direction. For example, on July 
20 the boat showed light winds that veered dramatically from the 
NE to SW at 1 PM local time. The WRF runs at both resolutions 
indicated a similar wind shift but closer to 4 PM local time. This will 
shift the ozone maximum by only a few hours. 

Although the 1.33 km resolution CMAQ model simulation 
closely matched the median for the 10-day period, on certain days 
the model was unable to predict ozone accurately, with a high bias 
shown in Fig. 7a and a low bias shown in Fig. 7b. During an ex- 
ceedance day (July 12), the model had a consistent 10-15 ppbv high 



Fig. 5. Median hourly ozone concentrations (ppbv) at the SRVx’s location and the 
closest upwind ground monitoring station from July 11, 2011 through July 20, 2011 as a 
function of time. 



Fig. 6. Median ozone concentrations (ppbv) at the SRVx’s location and at the closest 
CMAQ. (1.33 km) grid point for each hour from July 11, 2011 through July 20, 2011 as a 
function of time. 


bias and on an earlier exceedance day (July 13), the model had a 
10-15 ppbv low bias. The high bias of the model can likely be 
attributed to the boundary layer depths calculated by WRF and 
input into CMAQ, while the low bias of the model may be related to 
a lower temperature at the surface or perhaps a more stratified PBL 
inhibiting downward mixing. 

Measurements of the aerosol-based boundary layer height were 
determined by a High Spectral Resolution Lidar (FISRL) instrument 
onboard the UC-12B aircraft on July 20. The HSRL dataset includes 
aerosol extinction at 532 nm, aerosol backscatter at 532 nm and 
1064 nm and depolarization at 532 nm and 1064 nm (Hair et al., 
2008) and profiles of aerosol backscatter are used to derive the 
mixed layer height (Scarino et al., 2013 ). Observations from the HSRL 
were compared to the modeled boundary layer from WRF on July 20. 
Observations were used only for July 20 because this was the only 
day the SRVx was in the north part of the Chesapeake Bay and the 
UC-12B aircraft simultaneously conducted a flight. On July 20, the 
modeled boundary layer in the morning agreed to within 100 m, but 
in the afternoon the modeled boundary layer was 300-500 m lower 
over the Chesapeake Bay than the observed aerosol-based boundary 
layer (Fig. 8). 

3.3. Observational comparisons: total reactive nitrogen 

Observations of reactive nitrogen species are critical since the 
eastern United States lies in the NO x -limited regime of ozone 

Table 3 

CMAQ model mean bias (model minus observations), normalized mean bias, root- 
mean square error (RMSE), and normalized mean error (NME) of ozone (ppbv) for 
the boat and nearby ground stations. 


Location 

Mean bias (ppbv) 

Normalized 
mean bias (%) 

RMSE (ppbv) 

NME (%) 

SRVx boat 

-0.78 

-1.55 

10.14 

20.23 

Southern MD 

5.00 

9.98 

10.69 

21.33 

Calvert County 

4.89 

9.77 

10.06 

20.08 

Seaford, DE 

2.76 

5.51 

9.64 

19.24 

PG Equestrian 

0.79 

1.57 

11.87 

23.68 

Davidsonville 

2.24 

4.47 

12.21 

24.35 

Beltsville 

1.60 

3.20 

10.52 

20.99 

Essex 

1.67 

3.33 

9.84 

19.63 

Edgewood 

-1.96 

-3.91 

11.05 

22.05 

Millington 

0.22 

0.44 

9.94 

19.84 

Lums Pond, DE 

1.77 

3.53 

12.40 

24.75 
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Table 4 

WRF model mean bias, max difference, and root-mean square error of 2 m tem- 
perature (°C) for the boat and BWI airport. 


Location 

Mean bias (°C) 

Max difference (°C) 

RMSE (°C) 

SRVx boat 

-0.52 

4.62 

1.59 

BWI airport 

1.24 

4.20 

2.35 


production (Chameides et al., 1992; Trainer et al., 1993; Frost et al., 

2006) due to the excess of largely biogenic isoprene. Accurate 
model output of NO y species is especially important due to reactive 
nitrogen’s critical role in ozone formation in the NO x -limited 
regime found in eastern United States during the summer. 

Observations from the SRVx were compared to the UC-Berkeley 
thermal dissociation laser-induced fluorescence (TD-LIF) instru- 
ment (Day et al., 2002) used on the P3-B when it flew spirals over 
the Chesapeake Bay. The TD-LIF does not measure NO, so all com- 
parisons are NO y — NO. The observations of NO y — NO from the 
SRVx using a chemiluminescence instrument with external mo- 
lybdenum converter are higher than the data from the TD-LIF. This 
is an expected outcome since NO y concentrations decrease expo- 
nentially with height (Brent et al., 2013) due to emissions that come 
from the surface and relatively short lifetimes compared to other 
trace gases. There were no other suitable ground observations of 
NO y during this campaign. 

3.4. Model comparisons: total reactive nitrogen 

Observations of NO y from the SRVx were compared to 1.33 km 
CMAQ results over the Bay. On each day of the 10 daytime cruises, 
with the exception of July 19 when the instrument was taken off- 
line for calibration, NO y observations were consistently lower 
than the output from the nearest grid point in CMAQ. The model 
regularly overestimated NO y and on July 12, it was overestimated 
by 100% in the mid-afternoon as shown in Fig. 9. The data from the 
TD-LIF instrument (Day et al., 2002) on the P3-B aircraft during a 
spiral on July 20 also indicate a significant overestimation of NO y 
species by CMAQ as shown in Fig. 10. While the vertical profiles of 
NO 2 and FINO 3 match well, alkyl nitrates (ANs) and peroxy nitrates 
(PNs) are overestimated by factors of 2 and 4 respectively. This 
overestimation of reactive nitrogen species has also been seen in 
other modeling studies (Brioude et al., 2013; Yu et al., 2012). 

To understand whether the overestimate is an emissions issues, 
chemistry issue, or both, we examined the partitioning of the NO y 



Fig. 7. Ozone concentration (ppbv) as a function of time at the SRVx’s location and at 
the closest CMAQ. grid point on a) July 12, 2011 and b) July 13, 2011. 


species. If partitioning is correct, then the issue is likely due to high 
emissions or low dispersion rates. To gain insight on this issue, we 
took the ratio of NO/NO y during the morning hours when the two 
species are positively correlated and the NO measurement is above 
the detection limit. As seen in Fig. 11, the NO/NO y ratios between 
the model simulation and observations often lie below the 1 -to-l 
line. The mean of the data shows NO concentrations are 10.0% of 
total NO y in the observations, while NO concentrations are 7.6% of 
total NO y in the CMAQ simulation. This indicates that CB05, as 
employed, partitions more NO y species as higher oxides (i.e., ANs, 
PNs, HNO3) than is observed. This suggests that gas-phase chem- 
istry scheme (CB05) overestimates the lifetimes of higher order NO y 
species such as ANs and PNs, deposition rates are too slow, or 
conversion rates of NO y to N0 2 are slower than observed. 

To minimize computing time, the CB05 chemical mechanism 
simplifies the alkyl nitrates by grouping all alkyl nitrates in a single 
chemical species (NTR). The lifetime of NTR calculated during a 
simulation of CMAQ using 2007 summer conditions, yields a life- 
time of 10 days. It has been shown that isopropyl nitrate has a 
lifetime of 10 days (Luke et al., 1989), but higher-order alkyl nitrates 
have a much shorter lifetime (1—2 days) (Horowitz et al., 2007; 
Perring et al., 2009), due to a lack of electronegativity holding the 
gas phase species together. The shorter lifetimes of the high-order 
alkyl nitrates species are not accounted for in the CB05 gas-phase 
chemistry scheme. After decomposition, the alkyl nitrates split 
into an alkyl chain and N0 2 . If the lifetime of NTR in CB05 were to 
be shorter, then this would yield lower concentrations of alkyl ni- 
trates, which would be more consistent with observations. 

To represent peroxy nitrates in the model, the CB05 mechanism 
simplifies the species into peroxyacyl nitrates (PAN), all other 
higher order peroxyacyl nitrates (PANX) and peroxynitric acid 
(PNA), with the latter being a very small fraction of the first two at 
high temperatures. The sum of the concentrations of peroxy ni- 
trates (PNs) in the model is higher than observed. The primary 
destruction of peroxyacyl nitrates is via thermal dissociation. At 
higher temperatures, PAN and PANX dissociate more rapidly into 
acetylperoxy radicals (CH 3 C( 0 ) 0 2 ) and higher order acylperoxy 
radicals (C 2 H 5 C(0)0 2 ) respectively. The concentration of PAN and 
PANX is therefore governed by the kinetic equilibrium rate con- 
stant, which is a function of temperature, and the concentrations of 
the products, CH 3 C(0)0 2 , C 2 H 5 C(0)0 2 , and N0 2 . There are stark 
differences in the kinetic equilibrium rate constant (JCeq) between 
IUPAC (2010) and JPL (2011 ), with the latter being 24% less than the 
former (3.03 x 10 -8 vs. 2.3 x 10 -8 ) at 298 K. The CB05 mechanism 
uses the higher IUPAC (2010) kinetic equilibrium rate constant, 
which favors a higher production rate of PAN. Furthermore, some 
studies (Turnipseed et al., 2006; Wu et al., 2012) have suggested 
that the dry deposition rates of PAN in the air quality models are too 
slow. Updating the rate constants of PAN formation as well as 
changing the dry deposition velocities, may better align the model 
output with observations. 

4. Discussion 

The observations from the SRVx show, with a 95% confidence 
level certainty, that ozone concentrations are elevated over the Bay 
when compared to upwind ground sites. The extended period of 
high ozone causes a larger number of days to exceed the U.S. EPA 8 - 
h 75 ppbv NAAQS threshold over the Bay than over nearby land 
areas. Here we discuss potential reasons for this phenomenon and 
attempt to apportion a relative importance for each mechanism. 

During a day that lacks precipitation, which is the case for most 
ozone exceedance days in the Baltimore-Washington metropolitan 
region, ozone is primarily destroyed by the mechanisms listed in 
Table 5 (Seinfeld and Pandis, 2006). Dry deposition (mechanism 1) 
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Fig. 8. PBL depth output by WRF minus measurements of boundary layer height using a high spectral resolution lidar (HSRL) aboard the UC-12 aircraft on July 20, 2011. 


is the primary mode of ozone destruction near the surface. Titration 
due to NO (mechanism 2) also occurs near the surface, but this 
serves as a reservoir for O 3 as N0 2 is re-generated. Mechanisms 3-5 
are most prominent in the upper troposphere and isolated ocean 
regions where dry deposition rates are minimal. 

Deposition is the primary mode of destruction in the boundary 
layer and occurs fastest in heavily forested areas (Fowler et al., 
2001; Nowak et al., 2006). Differences in ozone dry deposition 
rates have been widely studied. A list of 24-h-averaged dry 



Fig. 9. Total NO y concentration measured on the SRVx compared to total NO y from the 
closest grid point in CMAQon July 13, 2011 as a function of time. 


deposition velocities from the literature is given in Table 6 (Wesely 
and Hicks, 2000; Chang et al., 2004; Nowak et al., 2006). 

For a mix of 50% deciduous forest, 25% grass, and 25% pavement, 
the 24-h averaged dry deposition velocity for ozone is 0.50 cm s' 1 . 
However, estimates for dry deposition of ozone in coastal envi- 
ronments are 0.15 cm s _1 . The slower deposition velocity is due to a 
lack of vegetation and surface roughness in coastal areas (Gallagher 
et al., 2001 ). To calculate the difference in ozone deposition over an 
hour, we can use the formula described in Table 7. For a boundary 
layer depth of 800 m, which is typical over the Chesapeake Bay 
during the mid-afternoon, ozone concentration would be approx- 
imately 1 . 6 % higher after an hour than an air parcel of similar 
concentration over land due to slower deposition velocities over 
water, assuming all other environmental conditions are the same. 

If winds are from the southwest, maximizing residence time 
over the Bay, an air parcel that entered the southwest portion of 
Chesapeake Bay may have been over the Bay for approximately 5 h. 
By the time an air parcel leaves the Bay, its ozone concentration 
theoretically could be 8 % higher than transport over land. 

Boundary layer height also plays a major role in determining 
concentrations of ozone near the surface (Rao et al., 2003). Pol- 
lutants are primarily confined within the boundary layer due to a 
strong subsidence inversion during anticyclonic events. The only 
mechanism by which pollutants can be vented out of the boundary 
layer during strong anticyclonic setups is through fair-weather 
cumulus clouds (Dacre et al., 2007). However, cumulus clouds are 
largely non-existent over the Chesapeake Bay during strong sub- 
sidence events (Loughner et al., 2011). 

The boundary layer over land tends to be deeper because the 
surface temperature is higher over land during clear-sky conditions 
in the mid-afternoon. As the boundary layer depth decreases, 
emissions of ozone precursors, such as NO x compounds, 
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Fig. 10. a. Total NO y concentration (minus NO) split by compound (N0 2 , peroxy ni- 
trates (PN), alkyl nitrates (AN), and HN0 3 measured on the P3-B as function of altitude 
during the 1630 UTC spiral on July 20, 2011 over the Chesapeake Bay. b. Same as 
Fig. 10a but using model results from the closest grid point in CMAQat 1700Z (1 PM 
local time). 


accumulate in a smaller volume of the atmosphere leading to 
higher concentrations. On July 20 between 20 and 21 UTC or 4 PM 
and 5 PM local time, the HSRL aboard the UC-12B aircraft measured 
the aerosol-based boundary layer depth to be 1000-1200 m over 
land and 400-600 m over the Chesapeake Bay within 10 min as 



Ratios from SRVx Observations 

Fig. 11. NO/NO y ratios from 1.33 km CMAQ. run vs. observations from the ship during 
the morning hours when NO and NO y are positively correlated and NO is above the 
instrument’s detection limit. 


Table 5 

Loss mechanisms for ozone in the lower troposphere. 


Mechanism 1 

Dry deposition 

Mechanism 2 

O 3 + NO — > 0 2 + N0 2 

Mechanism 3 

0 3 + OH -*• 0 2 + H0 2 

Mechanism 4 

0 3 + H0 2 -► OH + 20 2 

Mechanism 5 

0 3 + hv -*■ 0 2 + 0( 1 D) at hv < 320 nm 


seen in Fig. 12. If there were no boundary layer venting and envi- 
ronmental conditions and emissions were identical, the concen- 
trations of NO 2 could be up to a factor of 2 higher over the Bay than 
over land leading to a substantial increase in O3, since the mid- 
Atlantic region is in the NO x -limited regime. However, there is 
likely some vertical mixing and emissions are likely lower over the 
Bay. Although there were no direct measurements of NO 2 at the 
surface of the bay during this particular campaign, data on the P3-B 
shows that at 0.3 km, the lowest altitude of the flight spirals, NO 2 is 
higher by as much as 0.5 ppbv over water than land. Using ozone 
efficiency rates from the DISCOVER-AQ campaign, for every 1 ppbv 
increase in NO x , ozone production will increase by an average of 
8.26 ppbv with a 90% confidence interval of 4.93-19.4 ppbv (He 
et al., 2013); this is slightly higher than an urban study in Hous- 
ton, which showed an average ozone production efficiency of 5.9 
(Neuman et al., 2009). 

A bay-breeze circulation often develops over the Chesapeake 
Bay during the late spring and early summer (Ryan et al., 1998; 
Stauffer et al., 2012) impacting the coastal temperature structure 
and associated meteorological conditions. The bay-breeze yields a 
meso-high pressure directly over the Chesapeake Bay, and a meso- 
low pressure just inland from the Bay. This creates stagnation and 
clear skies directly over the Bay. Fewer cumulus clouds develop 
over the Chesapeake Bay than over land because of the lower sur- 
face temperature, shallower boundary layer depth and relative lack 
of thermals over the water. Decreased cloud cover increases 
photolysis rates by allowing more UV radiation to reach the lowest 
levels of the atmosphere creating an environment more favorable 
for ozone production. On July 20, visible satellite imagery, seen in 
Fig. 13, shows an expanse of low level fair cumulus clouds over the 
Baltimore— Washington region, with no clouds over the Bay. Cloud 
coverage is estimated to be 10-30% over land and 0% over the Bay 
leading to a higher j(N02) value over the Bay. 

During the DISCOVER-AQ campaign, the P3-B aircraft measured 
j(N02) throughout its flights. In the mid-afternoon, 3:30 PM local 
time, on July 20, 2011 when the P3-B flew at an altitude of 390 m 
over land in an absence of clouds, the j(N 02 ) rate constant was 
0.0082 s -1 , while 30 s later underneath a fair-weather cumulus 
cloud, which was confirmed by looking at the forward camera on 
the P3-B, the j(N02) rate constant dropped to 0.0043 s _1 . If we 
assume the sky over land is filled with 20% cumulus clouds and the 
sky over the Bay has no clouds, the average j(NC> 2 ) would be 
0.0074 s _1 over land and 0.0082 s' 1 over the Bay. Therefore, 
dissociation of NO 2 into NO and odd oxygen may be up to 10.5% 
faster during the mid-afternoon of a summer day. 

It is estimated that NO x emissions from barges that travel the 
Chesapeake Bay account for 10% of all mobile emission sources (U.S. 


Table 6 

Ozone deposition velocities for various surface types. 



Forest 

Coastal 

Ocean 

Nowak et al. (2006) 
Wesely and Hicks (2000) 
Chang et al. (2004) 
Gallagher et al. (2001) 

0.5 cm s -1 
0.8 cm s -1 

0.148 cm s ” 1 

0.05 cm s " 1 
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Table 7 

Calculation of the difference in ozone dry deposition over land and the bay. 
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yj(land) 

PBL depth 

18 m/h _ 
800 m 


^d(bay) 

PBL depth 

5.328 m/h _ ro u~ 
800 m - 1 -* 3 ° n 


EPA, 2010). In March 2010, the U.S. EPA adopted a regulation 
requiring large barges to burn cleaner fuel that which emits less 
NO x when they are within 200 nautical miles of the North American 
coastline (U.S. EPA, 2010). However, this regulation was not 
enforceable by the U.S. EPA until August 2012, which is after the 
Maryland DISCOVER-AQ field study. Many large transport tankers 
burn bottom-of-the-barrel bunker fuel, which releases a higher 
proportion of NO x than diesel fuel (Eyring et al., 2005). To date, 
there has been little quantification of barge emissions (Mason et al., 
2008). Using the 8.26 ppbv O 3 per ppbv NO x ozone production ef- 
ficiency calculated during the DISCOVER-AQ campaign (He et al., 
2013), we estimate that 0.1 ppbv increase in NO x concentrations 
over the Chesapeake Bay could yield a 0.8 ppbv increase in ozone, 
since the mid-Atlantic region is characterized by the NO x -limited 
regime of ozone production. 

Halogen chemistry may play a role in ozone formation over the 
Chesapeake Bay. Recent modeling studies suggest that CI 2 photo- 
chemistry may result in an increase of 5-8 ppbv in daily maximum 
ozone levels (Finley and Saltzman, 2006). To see if more chlorine is 
available over the Bay, we looked at the 5 -year average (between 
2007 and 2011 ) of Cl" dry and wet deposition at two Clean Air Status 
and Trends Network (CASTNET) sites in Maryland. The Blackwater 
National Refuge site is located on Maryland’s eastern shore and is 
generally downwind of the Bay, while the Beltsville site is located 
upwind of the Bay. Dry deposition rates of Cl - are 2.14 times higher 



Fig. 13. Visible image from the MODIS satellite at 1610Z (2:10 PM local time) on July 
20, 2011 showing the presence of low-level cumulus clouds only over the land. 


over a 5-year average at the Blackwater site and wet deposition 
rates of Cl - are 3.62 times higher. 

One factor inhibiting ozone production over the Bay is the lower 
tropospheric temperature profile. Coastal areas in extratropical 
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Fig. 12. Measurements of boundary layer height using a high spectral resolution lidar (HSRL) aboard the UC-12 aircraft on July 20, 2011. 
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latitudes heat up more slowly than nearby inland locations during 
the summer due to the influence of the cooler waters. During the 
10-day campaign, temperatures on the SRVx at 2 PM local time 
were on average 3.4 °C cooler than the Baltimore-Washington In- 
ternational (BWI) airport which is located 30 km inland from the 
Chesapeake Bay. 

The dissociation of PAN into NO 2 has a strong temperature 
dependence (Seinfeld and Pandis, 2006). A calculation of the rate 
constant using IUPAC (2010) shows that PAN dissociates 1.66 times 
quicker at 304.3 K than 300.9 K. The quicker dissociation of PAN at 
higher temperatures over land shifts the equilibrium reaction to- 
ward N0 2 , the primary precursor to ozone in the NO x -limited 
regime over the Mid-Atlantic. However, the dissociation of PAN is 
slower over the Bay, keeping more N0 2 tied up as PAN, and thereby 
decreasing O3 production. 

5. Conclusions 

Observations from the NOAA SRVx vessel during the DISCOVER- 
AQ and GEO-CAPE CBODAQ campaigns show with a certainty 
exceeding the 2-sigma level, that daytime ozone concentrations are 
elevated over the Bay when compared to the closest upwind 
ground station. We posit that this high anomaly is influenced by a 
number of mechanisms, in approximate descending order: 

• Shallower boundary layers trapping shipping emissions near the 
surface 

• Higher photolysis rates due to clear skies over the bay 

• Decreased boundary layer venting due to a lack of fair-weather 
cumulus clouds 

• Slower deposition velocity over the Bay 

The ozone concentrations exhibit a high anomaly over the Bay 
even though temperatures are cooler and allow precursors to ozone 
such as PAN to remain more stable. The observed high anomaly 
over the Chesapeake Bay is of primary importance since many 
citizens spend their leisure time on or near the Chesapeake Bay 
during the summertime, and are exposed to the unhealthy air 
quality conditions. Onshore winds can bring these pollutants to 
local coastal and inland communities. Expanded monitoring of 
ozone directly over the Chesapeake Bay is needed to precisely 
quantify the extent of this high anomaly. 
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